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Preface 

This  report  demonstrates  how  receiver  operating  characteristics 
may  be  calculated  for  a broad  class  of  nonlinear  detectors,  those 
whose  outputs  are  quadratic  forms  of  the  inputs.  Of  particular 
Interest  is  the  capability  to  perform  these  calculations  for 
arbitrary  bandwidth-integration  time  products  and  dependent 
time  samples.  The  work  documented  by  this  report  was 
performed  in  the  Signal  and  Digital  Processing  Branch  under 
Task  No.  A03S370B/001B/7F11/100-000. 


EDWARD  C.  WHITMAN 
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COMPUTING  R.O.C.  FOR  QUADRATIC  DETECTORS 
Prepared  by: 

L.  E.  Miller 

CHAPTER  1 

COMPUTING  RECEIVER  OPERATING  CHARACTERISTICS 
FOR  QUADRATIC  DETECTORS 

1.  INTRODUCTION 

For  comparing  performances  of  different  signal  detector  configura- 
tions, the  functions  known  as  receiver  operating  characteristics  (ROC) 
are  a standard  tool.  A given  detector's  ROC  may  be  expressed 

Pp  = Y(h2;o)  (1-1) 

where 

Y = Pp  = probability  of  detection 

a i Pp^  = probability  of  false  alarm 

h^  = input  slgnal-to-noise  ratio  (SNR).  (1-2) 

The  probabilities  are  computed  on  the  assumption  that  the  detector 
output  undergoes  a statistical  test  to  decide  between  the  hypotheses 

H^ : no  signal  is  present 
H^ : a signal  is  present 

The  test  is  performed  by  comparing  the  value  of  the  detector  output 
z to  a number  t called  a threshold.  If  z ^ t,  H^  is  accepted;  if 
z < T,  H is  accepted. 
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If  there  actually  Is  no  signal  (h  *0)  and  z ^ t,  then  the  test 
yields  an  error  of  Type  I or  false  alarm,  with  the  probability 

w 

o - y dz  Pj^(z|  h2»0)  (1-4) 

T 

where  P_(z|h^=0)  is  the  probability  density  function  (pdf)  of  the 
z 

detector  output  subject  to  the  condition  that  there  is  no  signal. 

If  the  value  of  t Is  chosen  to  keep  o at  a particular  value,  then  a 
is  known  as  the  "level"  of  the  test. 

If  there  actually  is  a signal  present  but  z < t,  the  test  yields 
an  error  of  Type  II,  with  the  probability 

6 = y dz  P2(zlhV0).  (1-5) 

Conversely,  if  z ^ 0 when  there  is  a signal  present,  the  test  result 
is  correct  with  the  probability 

CO 

Y = 1 - 0 = y dz  P2(z|h2?^0)  . (1-6) 

T 

If  T has  been  constrained  by  o so  that  t = T(a),  then  y = Y(h^;o)  and 
is  known  as  the  "power  function"  of  the  test. 

It  is  clear  that  to  compute  the  ROC  for  a given  detector  requires 
knowledge  of  the  probability  distribution  of  the  detector  output. 

Often  the  functional  form  of  this  distribution  is  very  difficult  to 
obtain  analytically. 

The  purpose  of  this  report  is  to  document  methods  for  calculating 
the  ROC  of  a class  of  detectors.  Attention  is  first  paid,  in 


6 


NSWC/WOL/TR  76-li>8 

Chapters  2-4,  to  cases  in  which  samples  of  the  detector  Inputs  are 
considered  Independent,  thus  allowing  direct  calculation  of  the  ROC. 
In  Chapters  5 and  6,  the  more  general  case  Is  treated,  using  an 
approximation  method. 

1.1  Channel  Model. 


It  Is  assumed  that  the  detector  has  one  or  more  Inputs  x. (t). 


where 


x^(t)  = s^(t)  + n^(t),  1=1,2, ...,L. 


(1-7) 


and  the  noise  terms  n^(t)  are  Jointly  Gaussian,  stationary  random 
processes  with  zero  means  and  LxL  covariance  matrix  with  elements 


E{n^(t)n^(t)}  = 


(1-8) 


Two  types  of  signals  will  be  considered,  deterministic  and  random. 
Random  signals  will  be  considered  to  be  from  Jointly  Gaussian, 
stationary  random  processes  with  zero  means  and  LxL  covariance 


matrix  with  elements 


E{s^2(t)}  = 


E{5^(t)Sj  (t) } = • 


(1-9) 


Bandwidth  and  spectra.  Two  categories  of  detector  bandwidth  will 
be  treated,  low-pass  and  narrowband.  In  either  case,  a (two-sided) 
bandwidth  of  W hertz  and  Ideal  (flat)  response  Is  assumed.  This 


F 
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bandwidth  Is  lumped  with  the  signal  and  noise  models,  so  that, 
assuming  flat  spectra  over  the  bandwidth,  we  have 


= Nj^W,  dj^2  = 


(1-10) 


for  the  random  waveforms.  Deterministic  signals  are  modeled  as 


Sj^Ct;  = S^(t)cos[a)^t-0^(t)]  (1-11) 

In  which  the  envelope  and  phase  are  slowly  varying,  so  that 
signal  power  Is  given  by 

s^2(t)  = s^2,  lowpass  signal 

Sj^2(t)  = S^2/2,  narrowband  signal.  (1-12) 

Thus  Input  SNR's  are  taken  to  be 

s^/o^,  lowpass  deterministic  signal 
S^/2o^t  narrowband  deterministic  signal 
d^/a^,  random  signal.  (1-13) 

An  alternate  form  for  the  Information  represented  In  a detector’s 
ROC  Is  to  write 

h2  = h2(Y,«),  (1-14) 

a function  which  answers  the  question,  "What  SNR  Is  required  to 
produce  Pp  “ Y when  the  false  alarm  probability  Is  a?"  Often  this 
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I 

J 


I relationship  Is  abbreviated  by  giving  its  value  for  specified 

t 

j (y,o)  - (yi  0i)»  say  (y.o)  = (.S.io"**),  and  calling  it  "minimum 

detectable  signal  (MDS): 

MDS  = h^CYijOii).  (1-15) 

Also,  for  deterministic  signals  a number  called  "detection  threshold" 
(DT)  is  often  quoted  and  is  the  MDS  referred  to  a one-hertz  bandwidth: 

DT  = Wh2(Yi,ai)  . (1-16) 

In  the  case  of  narrowband  signals  and  noise,  the  detector  Inputs 
can  be  written 

Xj^(t)  = X^(t)cos[u)Qt-((.j^(t)] 

= (t  )cosa)Qt  + x^g(t)slna)^t 

= (n^^^+s^^  )cosiD^t  + (n^g+Sj^g)sin(»)Qt.  (1-17) 

Under  this  expansion,  n^^  and  n^^  are  independent  and  each  set  ^’^ic^ 

and  {n,  } has  covariance  matrix  T . A similar  statement  is  true  of 
is  n 

the  signal  terms  if  they  are  Gaussian. 

The  envelopes  of  the  random  Inputs'  correlation  functions  have 
the  factors 

j sln(irWAt  )/irWAt , narrowband 

sln(2TrWAt  )/TrWAt , lowpass  • (l-l8) 
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Therefore,  if  samples  of  the  Inputs  are  taken  1/W  and  1/2W  seconds 
apart,  respectively,  the  samples  are  uncorrelated  and,  being  Gaussian, 
Independent . 

1.2  Detector  Model. 

Let  M samples  of  the  L detector  Inputs  be  assembled  to  form 
a vector: 


C'  ~ CXi(ti),...,Xy(ti)}Xi(t2),.».,  Xj  (t2)j...jX]^(tj^),...,  X^  ^ ^ * 


(1-19) 


where  the  prime  (')  Indicates  transpose.  Detectors  will  be  considered 
which  produce  as  an  output  decision  variable  or  statistic  the 
quadratic  form 

z = LMxLM,  (1-20) 

and  therefore  termed  "quadratic  detectors."  In  Part  II  (Chapters  5,6), 
expressions  will  be  developed  at  this  level  of  generality  In  order 
to  treat  correlated  Input  samples.  However,  In  Part  I (Chapters  2-4) 
and  In  calculations  based  on  the  more  general  cases,  It  will  assumed 
that 

M 

z = I ^5’jO.Cj,  Q LxL  (1-21) 

where  ^ ^ 

5’j  = Cxi(tj),  X2(tj),...,  Xj^(tj)].  (1-22) 
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Effectively,  this  Is  the  same  as  saying  that  Is  a block  diagonal 
matrix  with  Identical  LxL  diagonals  Q/M. 

The  detector  form  (1-21)  Is  Illustrated  In  Figure  1-1.  The 
"quadratic  processor"  Is  assumed  memoryless,  so  that  the  sampling 
operation  can  be  placed  after  the  processor  for  convenience. 

This  model  Is  of  Interest  because  It  Is  an  Idealized  representation 
of  the  non-sampled  case  with  post-detection  Integration. 

That  Is, 


1 

M 


M 

I z(jAt) 
J-1 


T / 


dt  z(t) 


(1-23) 


If  At  » 1/W,  then  the  equivalence  M ■ WT  Is  made,  and  Is 

commonly  referred  to  as  the  "tlme-bandwldth  product"  In  communications 

theory. 
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CHAPTER  2 


COHERENT  POWER  DETECTOR 


In  this  chapter,  and  in  the  two  which  follow,  three  types 


of  quadratic  detector  of  frequent  Interest  are  treated  for  the 


case  in  which  time  samples  {ZjJ  of  the  processor  output  are 


Independent,  allowing  us  to  compute  ROC  exactly. 


2.1  Detector  Configuration. 


A great  variety  of  practical  detectors  can  be  modeled 


by  the  configuration  of  Figure  2-1,  in  which  the  detector  Inputs 


x^(t)  are  summed  and  then  squared.  Commonly  it  is  called  the 
"square-law  detector",  and  sometimes  the  "conventional  detector." 


Let  the  sum  of  inputs  be 


x(tj)  = I x^(tj) 


(2-1) 


NSWC/WOL/TR  76-1^8 


80  that 


E{x} 


8 ■ ^2^1 

0 RS 


(2-2) 


and 


Var{x} 


DS 

j,2+d2  -o2+  RS 


(2-3) 


Thus,  for  example,  o^»Oj^+2pOja2+02^  ^ deterministic 

signal . 

In  the  lowpass  input  case,  the  detector  output  is 


M 


M ^ * 

J“1 


M-2WT. 


Now,  the  sum  of  the  squares  of  v independent  Gaussian  random  variables 

with  unit  variances  and  means  is  a noncentral  chi-squared  variable 

2 

(x'^)  with  V degrees  of  freedom  and  noncentrality  parameter  X » * 

Therefore,  we  have 


2WTz 


is  X'2{2WT,2WTh2}  DS 


and 


1^  is  x2{2WT}  , 


RS 


(2-5) 


(2-4) 


where  x^  denotes  a (central)  chi-squared  variable  - one  for  which  x«0. 
•Using  "DS"  for  deterministic  signal  and  ”RS"  for  random  signal. 
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In  the  narrowband  input  case  the  detector  output  Is 

M 

^ ^ ^ m-wt 

J“1 


(2-6) 


Since  ■ S^^cose^  and  E{x^^}  * S^sine^  for  deterministic  signals, 

we  have 


is  x'^(2WT,2WTh2}  DS 


2WTz  . 2,^,™, 

is  X {2WT}  , 


the  same  result  as  for  the  lowpass  case. 
2.2  Receiver  Operating  Characteristics 
From  (1-6)  we  have 


RS 


(2-7) 


■/ 

■/ 


dzp^(z) 


dv  p^,2(vl2WT,2WTh2) 

2^r/a^ 

= Q(2WTT/a2 |2WT,2WTh^),  DS 


(2-8) 


where  Q(x*^|v,X)  is  the  noncentral  chi-squared  probability  Integral. 
Similarly,  we  have 


Y = Q[2WTT/(a2+d2) |2WT],  RS 

where  Q(x^lv),  the  chi-squared  probability  Integral,  is  widely 
tabulated. 


(2-9) 


16 


NSWC/WOL/TR  76-Uii 

EXAMPLE  2-1.  Problem:  What  Is  the  probability  that  a chl-squared 
random  variable  with  18  degrees  of  freedom  exceeds  the  value  30? 
Solution:  Using  Table  26.7  of  [1],  we  find  that  Q(30|l8)  = .03745. 
Setting  h^-o  In  (2-8)  or  d^-O  In  (2-9)  yields 


P = o = Q(2WTT/a2|2WT). 


(2-10) 


EXAMPLE  2-2.  Problem:  Find  the  threshold  value  t necessary  to  maln- 

-4 

tain  a false  alarm  probability  of  10  when  WT  = 100. 

Solution:  In  Table  II  of  [2],  the  function  P(a,c)=l-Q(2a/2c ) Is 

tabulated.  Under  P=.9999  and  c=100  we  find  that  a=l4l.530  or 

T«2a2a/2WT*1.4l53o2. 

The  procedure  for  obtaining  ROC  Is  diagrammed  In  Figure  2-2. 

For  either  type  of  signal,  the  process  begins  with  the  selection 
of  a false  alarm  probability  a.  Using  the  Identity 


and  (2-10),  the  false  alarm  threshold  Is  found  to  be 


■^a  “ 2WT  Xj_o(2WT)  . 


(2-11) 


(2-12) 


We  may  speak  also  of  a normalized  false  alarm  threshold  d^,  given  by 


- E{z}  Xi_o(2WT)-2WT 


(2-13) 


2/WT~ 


ONE  POINT  ON  ROC 

FIG.  2-2.  PROCEDURE  FOR  COMPUTING  ROC 
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since  E{x^(v)}  *>  v and  Var{x^(v)}  *2v.  The  quantity  Is  Interpreted  ’ 
as  the  number  of  detector  output  noise  standard  deviations  above  the  | 

mean  at  which  the  threshold  t must  be  set  to  produce  a false  alarm  j 

probability  a. 

EXAMPLE  2-3.  Problem:  Find  d for  WT»10  and  a-. 01. 

— ■ ^ a 

Solution:  The  table  on  page  252  of  [3]  gives  P=Q(x^|v).  For  P«.01 
and  v-20,  we  find  x^qo(20)-37.566.  Therefore  d =(37.566-20)/2/TD-2.777. 

For  convenience,  values  of  x^_j^(2WT)  and  d^  are  given  In  Table  2-1 
for  various  values  of  a and  WT.  Also  d^  Is  plotted  vs  a for  WT=1,  10, 

50,  -In  Figure  4 of  [4]. 

Having  fixed  the  value  of  t^,  the  ROC  Is  a plot  of  Pp=v(h^)  or. 

If  It  Is  more  convenient,  h^-h^Cy).  For  Gaussian  signals,  from  (2-9) 
and  (2-12)  we  have 

Y = QCx?  (2WT)/(l+h2)  |2WT].  RS  (2-14) 

1— a 

One  approach  Is  to  Interpolate  In  a chl-squared  table  to  find  the  value 
of  Q(x^|v)  corresponding  to  the  number 

X?  „(2WT) 

x2  (2WT)  = -i-2 . (2-15) 

1 ^ l+h2 

However,  It  Is  simpler  to  calculate  Q,  since  It  Is  given  by  the  finite 
sum  [1] 

Q(x2|2v)  = ''^\xV2)Vk! 

k=0 

= e-X^/2  e^_^(x2/2)  , (2-16) 
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/ 

/ 


i 


wt\ 

10'^ 

10'2 

10® 

lO’^ 

10-® 

A'va  <2WT)^ 

1 

4.60517 

9.21034 

13.816 

18.421 

27.64 

2 

7.77944 

13.2767 

18.467 

25.513 

33.38 

5 

15.9871 

23.2093 

29.588 

35.564 

46.86 

10 

28.4120 

37.5662 

45.315 

52.386 

65.42 

20 

51.8050 

63.6907 

73.402 

82.062 

97.66 

50 

118.498 

135.807 

149.449 

161.319 

182.12 

100 

226.0210 

249.446 

267.458 

283.060 

309.84 

do,  NORMALIZED  FALSE  ALARM  THRESHOLD^^ 

1 

1.3026 

3.6052 

5.9080 

8.2105 

12.82 

2 

1.3362 

3.2798 

5.1149 

6.8989 

10.3874 

5 

1.3388 

2.9537 

4.3800 

5.7163 

8.2421 

10 

1.3301 

2.7775 

4.0027 

5.1207 

7.1815 

20 

1.3198 

2.6487 

3.7345 

-.7027 

6.4466 

50 

1.3080 

2.5319 

3.4966 

4.3359 

5.8068 

100 

1.3011 

2.4723 

3.3729 

4.1530 

5.4920 

OO 

1.28155 

2.32635 

3.09023 

3.71902 

4.7534 

NOTES: 

•VALUES  FOR  WT  = 100  AND/OR  a = 10  ® TAKEN  FROM  TABLE  II  OF  [2]  ; 
THE  REMAINDER  FROM  TABLE  26.8  OF  [ 1 ] . 

••COMPUTED  USING  VALUES  OFX^va^^WT)  FROM  TABLE  II  OF  [2] . 

TABLE  2-1.  FALSE  ALARM  THRESHOLDS. 
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using  the  program  of  Figure  2-3.  This  calculation  was  performed  for 

a number  of  cases,  and  the  results  are  plotted  In  Figures  2-4  and 
_5  -4 

2-5  for  a^lO  and  10  , respectively. 

EXAMPLE  2-4.  Problem:  - What  Integration  time  T Is  required  to 
achieve  90%  detection  ahd  1%  false  alarm  probabilities  when  the  Input 
SNR  Is  zero  dB,  If  the  bandwidth  Is  0.25  Hz? 

Solution:  Prom  Figure  2-4,  y“.9  and  h^-O  dB  correspond  approximately 
to  WT=25.  Thus  about  100  seconds  of  Integration  time  Is  required. 

For  deterministic  signals,  from  (2-8)  and  (2-12)  we  have 


I 


1 


Y = Q[x^  (2WT) |2WT,2WTh2].  (2-17) 

l-a 

Here,  direct  calculation  Is  not  as  simple  since  the  noncentral  chl- 
squared  probability  Integral  Is  given  by  the  Infinite  sum  [1] 

Q(x'21v,X)  = e"^/^  I Q(x*^|v+2k).  (2-18) 

k=0  k! 

Tables  for  this  function  are  relatively  rare  [5,6],  of  which  [6]  Is 
fairly  extensive  for  2WT  ^ 100  and  o >,  10~^.  For  WT=1,  Q Is  well 
known  to  communications  engineers  as  Marcum’s  Q-functlon,  thoroughly 
tabulated  In  [7]  and  given  by  nomograph  In  [8].  Urkowltz  [9]  uses  a 
nomograph  for  the  (central)  chl-squared  probability  Integral  [10] 
and  an  approximation  technique  due  to  Patnalk  [11]. 

In  order  to  avoid  both  the  labor  of  Interpolating  tables  and  the 
Inherent  Inaccuracy  of  nomographs,  an  excellent  Gaussian  approximation 
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76/09/10.  14.22.33. 

PROGRAM  CHI534 


10  READ  A,  X,  M 

20  PRINT  "M  = "M,"A  ="A,"X  = "X 
30  PRINT 

40  PRINT  ”SNR(DB)",’'PD" 

50  FOR  V=2.5  TO  3 STEP  .1 
60  H=10+V 
70  Y=.5*  X/(l+H) 

80  S1=S2=1 
90  FOR  K=1  TO  M-1 
100  S2=S2*Y/K 
110  S1=S1+S2 
120  NEXT  K 
130  S1=S1»EXP(-Y) 

140  PRINT  10*V,S1 
150  NEXT  V 
160  STOP 
170  DATA  .0001 

175  data  23.5:3,2 

180  END 


H = h2 


SI  = Q(x2_y|2M) 

FIGURE  2-3 

Program  to  Compute  the  Chl-Squared  Integral 
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I 


■10  0 10  20  30 

INPUT  SIGNAL-TO-NOISE  RATIO  (h^),  dB 

FIG.  24  RECEIVER  OPERATING  CHARACTERISTICS,  SQUARE-LAW  DETECTOR, 
FOR  PFA  > 0.01  AND  WT  VARIED  (RANDOM  SIGNAL) 


I 


I 

i 


I 
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INPUT  SIGNAL-TO-NOISE  RATIO  (h^),  dB 

FIG.  2-5  RECEIVER  OPERATING  CHARACTERISTICS,  SQUARE-LAW  DETECTOR, 
FOR  PFA  > 0.0001  AND  WT  VARIED  (RANDOM  SIGNAL) 
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to  the  noncentral  chl-squared  distribution,  due  to  Sankaran  [12], 
was  used  to  obtain  in  terms  of  y In  Appendix  A it  is  shown  that 


h2(Y) 


A.  , ^ , 2wt^i 

/vrr  aitV  /w 


(2-19) 


where  d is  defined  by 

y 

y = Q(d^)  - |[l-erf(d^//D]  , (2-20) 

Q(x)  being  the  Gaussian  probability  integral*  and  erf(x)  the  error 
function,  both  universally  known  and  tabulated.  For  reference. 

Table  2-2  provides  several  values  of  d^. 

It  is  interesting  to  note  that  if  the  value  of  h^  yielding  y=.5 
is  taken  to  be  the  minimum  detectable  signal  (MDS),  then  from  (2-19) 
it  is  given  simply  as 


MDS  = 


(2-21) 


This  is  to  be  compared  with  the  heuristic  value  of  MDS=d^//WT  in  [4], 
which  is  based  o’^  distributional  properties  when  WT  is  very  large 
(l.e.,  the  central  limit  theorem). 

EXAMPLE  2-5.  Problem:  Find  the  MDS  for  a»10"^  and  WT>=50. 

Solution:  Since  no  table  of  y is  available  for  this  value  of  o, 

we  use  the  approximation  (2-21),  yielding  MDS“4 . 3359/»^5o+.  01=.  6232= 

-2.05  dB,  where  d was  obtained  from  Table  2-1.  If  W=.l  Hz,  the 

•One  of  the  spinof^fs  from  this  study  is  the  discovery  that  Marcum’s 
Q-functlon  is  rather  easily  and  accurately  approximated  by 

Q(a,B)  « Q(/3-l/2  - /o+l/2). 
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i 

I 


t 


r 

dy 

10-® 

4.75342 

10-® 

4.26489 

3.71902 

.0002 

3.54008 

.0005 

3.29053 

.001 

3.09023 

.002 

2.87816 

.005 

2.57583 

.01 

2.32635 

.02 

2.05375 

.05 

1.64485 

.1 

1.28155 

.2 

0.84162 

.3 

0.52440 

.4 

0.25335 

.5 

y-Q(dy)  I 

NOTES: 

(1) d,.y--dy  I 

(2)  VALUES  TAKEN  FROM  TABLES  26.5  | 

AND  26.6  OF  [1] . 


TABLE  2-2.  GAUSSIAN  PROBABILITY  INTEGRAL. 
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corresponding  detection  threshold  Is 

DT  - (.1)  MDS  - -12.05  dB. 

Using  (2-19),  ROC  were  computed  for  the  square-law  detector 
with  deterministic  input  signal,  and  are  displayed  in  Figures  2-6 
and  2-7  for  a ■ 10”^  and  10“**,  respectively. 
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PROBABILITY  OF  DETECTION  ( 7 ) 
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FIG.  2-7.  RECEIVER  OPERATING  CHARACTERISTICS,  SQUARE-LAW  DETECTOR, 
FOR  PFA  - 0.0001  AND  WT  VARIED  (DETERMINISTIC  SIGNAL) 
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I 


29 


NSWC/WOL/TR  76-148 


r 

I 


CHAPTER  3 

CORRELATION  DETECTOR 


3.1  Detector  Configuration 

A second,  broad  class  of  practical  detectors  can  be  modeled 
by  the  configuration  of  Figure  3-l»  In  which  the  detector  inputs 

1»1,2,...L}  are  first  partitioned,  then  the  partitions  summed 
and  multiplied  together.  Thus  the  output  y(t)  of  the  multiplier 
has  the  form 


’ LI 

L 

y(t)  ■ u (t)u  (t)  » 

1 2 

I 3C4(t) 
1-1  ^ 

1 Xi(t) 

1-Ll+l  ^ 

(3-1) 


Let  us  define  for  J*l,2 


E{ 


Uj(t))  . I 


Sj(t) 


DS 

RS 


(3-2) 


Var{Uj (t) } 


DS 
2 RS 


(3-3) 


and 


Cov{u  (t),u  (t)} 
1 2 


(P®1®2 

(po^o^+rd^d^. 


DS 

RS 


(3-4) 


30 
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The  product  can  be  written 


u u 
1 2 


u2  - u2 
3 *« 


where 


u “ i(u  + u ) 


u » tCu  - u ) 
4 1 2 


The  two  new  Gaussian  variables  thus  defined  have  the 


E{u  } 
3 


S “ t(s  +s  ) 
3 12 


E{u  } 
4 


S ■ ■f(Si-S2) 


DS 

RS 

DS 

RS 


Var{u  } 
3 


o2  =-J(a2+2po  0 +o2) 
3 1 12  2 


o2+d2-o2+^(d2-2rd  d +d2) 
3 3 3 1 1 2 2 


Var{u  } 

4 


r2®i(o2+2po  o +0^) 


1 2 2 


o2+d2*a2+r (d2— 2rd  d +d2) 

4 4 4 1 1 2 2 


DS 

RS 

DS 

RS 


and 


Covfu  ,u  } 
3 4 


/no  0 *5(o2-o2) 
3 4 12 


no  o +6d  d «i(o2+d2-o2_d2 j 
3 4 3 4 ^1122 


In  the  lowpass  case  the  detector  output  Is 


32 


(3-5) 


(3-6) 

moments 


(3-7) 


(3-b) 


DS 

RS.  (3-9) 
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M 


^ I Cu2(ti)  - U2(t.)] 


A z - 2 , M-2WT 

“ 1 2 


where  for  deterministic  signals 


2WTz 


is  x'^(2WT,2WTh2) 

|2  3 


DS 


2WTz 


is  x'^(2WT,2WTh2) 

i2  4 


and  for  random  signals 


2WTz 


o*  + d^ 
3 3 


and 


2WTz 


2 j.  /i2 

o + d 
4 4 


are  x (2WT).  RS 


(3-10) 


(3-11) 


(3-12) 


Thus  the  detector  output  is  the  difference  of  two  scaled,  correlated 

(central  or  noncentral)  chi-squared  random  variables.  However,  from 

(3—9),  if  the  input  noise  power  is  the  same  in  each  channel  (o  ®o), 

1 2 

these  output  variables  are  uncorrelated  (and  Independent)  for  determlnls- 

2 2 

tic  signals;  for  Gaussian  signals,  the  same  effect  results  if  o + d » 

0*  + d^. 

2 2 

In  the  narrowband  case  the  detector  output  is 

33 


NSWC/WOL/TR  76-148 


2 1,22 

h * [h  +h  +2h  h cos(e  -e  )] 

3 2U+P;  12  11  12 

DS 

2 T 2 2 

h - V//  ^[h  +h  -2h  h cos(0  -0  )].  (3-14) 

4 2(lnp  ) 1 2 12  12 

3.2  Probability  Distribution. 

The  functional  form  of  the  distribution  of  the  correlation 
detector  output  z.  In  the  case  of  no  post-detection  Integration  (M“l), 
has  been  developed  In  complete  generality  for  the  assumed  model  [13,19]. 


This  form  Is  very  complex  and  far  from  being  "closed".  A great 

i 
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simplification,  however,  occurs  if  the  equal  input  variance  conditions 

discussed  above  are  assumed,  so  that  z and  z are  independent. 

1 2 

Using  those  assumptions,  in  Appendix  B it  is  shown  that  the 
probability  density  function  of  z becomes,  for  deterministic  signals, 

2 M-1  2 2 Y 

p(z)  - JL(i^£_)  exp  |-M(h  +h  )} 

34; 

M 00 

I I ;rCTCM(l-p)hV2j”CM(l+p)hV2]'"  DS  (3-15) 

(M=WT) 

exp  j-2Mz/a^(l-p)  j GjJ;[;”“J[4Mz/aNl-p^)J,  z>  0 

X 

exp  j2Mz/a^(l-p)|  GjJJ;J}:];[-4Mz/a^  (1-p^  ] , z<0 

where 

G^'Cx)  t I (3-16) 

n k=o  ® 


1 


are  polynomials  whose  properties  are  discussed  in  [13]*  For  Gaussian 

2 2 

signals,  the  pdf  is  obtained  by  setting  h =h  =0  in  (3-15).  Using 

34 


1 

! 


where  the  incomplete  exponential  function  Is  given  in  (2-16). 

Thus  the  probability  of  detection  for  the  correlation  detector  is 

,,  , 2 M-1  22  2 

Pd(-')  = exp{-M(h^+h^)-2MT/o  (1+p)} 


DS 

^ ^ =n[^CM(l-p)hy2]^[M(l+p)hV2]"‘ 

n-0  m*0  3 4 

’^‘'■5?“^  ,2M+n+m-2-K,i  , 2 J 

I m+M-1  ^ 4-p^  e„[2MT/a2(l+p)]  (3-20) 

K=0  ^ 

(M  - WT) 


36 


NSWC/WOL/TR  76-148 


! 

f 

j 

3“3  Receiver  Operating  Characteristics. 

* The  computer  program  presented  in  Figure  3-2  was  used  to 

calculate  a truncated  version  of  (3^20).  False  alarm  thresholds  are 

first  computed  according  to  the  procedure  diagrammed  in  Figure  3-3* 

-4 

which  Illustrates  the  case  of  o*10  and  p»0.  A number  of  values 
of  the  threshold  thus  obtained  are  listed  in  Table  3-1  for 
different  o and  WT. 

An  impression  of  the  rate  of  converge  of  the  distribution  of  the 
detector  output  to  that  of  a Gaussian  variable  can  be  gained  from 
Table  3-2,  in  which  thresholds  from  Table  3-1,  normalized  by 

✓VarCz}  = a^/(l+p)/2M  , (3-22) 

are  compared  to  their  Gaussian  counterparts  taken  from  Table  26.6 
of  [1]. 

Having  computed  false  alarm  thresholds  corresponding  to  various 
values  of  a,  the  method  of  Figure  3-4  was  employed  to  compute  Pp=Y 
for  fixed  a and  WT,  and  for  input  SNR  varied  when  the  signal  is 
deterministic.  For  computation,  the  case  for  which  there  is  no 
inter-channel  noise  correlation  (p=0)  and  for  which  the  signals  are 
of  equal  phase  and  power  (S  =S  =S,  0 =6  ) was  selected.  For  this 

O 12  12 

2 2 2 2 

case,  h =2h  and  h =0,  where  h is  the  input  SNR  in  each  channel. 

3 «♦  , 

— 2 

The  resulting  ROC  are  given  in  Figures  3-5  and  3-6  for  a=10  and  10  , 

respectively. 

EXAMPLE  3-1 » Problem;  Find  the  input  SNR  required  to  achieve  a 
letectlon  probability  of  .95  when  a is  constrained  to  be  1%  for  a 
correlation  detector  whose  bandwidth  is  0.1  Hz  and  the  integration 
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H3=2»H/(l+RO) 

H4  = 0 

IF  H3>0  THEN  73 
L2  = 0 

IF  H4>0  then  75 

LI  = 0 

C9=CO 

U0=.5*M*(1+R0)»H4 

VO=.5*M*(l-RO)*H3 

Y=2»M*X/(l+RO) 

S=0 

Ul=l 

CO=C9 

FOR  N=0  TO  LI 
IF  N=0  THEN  180 
U1=U1»U0/N 

CO=CO* ( 2«M+N-2 ) /N+M-1 ) 

Vl=l 

CL=CO 

FOR  J=0  TO  L2 
IF  J=0  THEN  240 
V1=V1*V0/J 

C1=C1* ( 2»M+N+J-2 )/ (M+J-1 ) 

W1=W8=1 

C2=C1 

FOR  K=0  TO  M+J-1 
IF  K=0  THEN  310 
W8=W8»Y/K 
W1=W1+W8 

C2=C2*(M+J-K)/(2«M+N+J-1-K) 

S=S+W1»U1*V1«C2*(2/(1-R0))+K 

NEXT  K 

NEXT  J 

NEXT  N 

S=S»EXP(-Y-M»(H3+H4)) 

S=S» ( 1+RO) * ( l-RO+2 ) + (M-1) 
S=S*2/4+M 


Mnemonics ; 


.Y  5 gi  , H3  = h2  HU  s hj. 


e.(Y),  Ul=[M(l+p)h2/2]^/N:, 

rMM  r?  - i2M+N+J-2-K 

tM(l-p)h^/2]  /J!,  C2  = \ ) 


(1+RO)  (1-R02)^‘^  exp  t-Y-M(H3+H4)] 


LI  L2  M+J-1 

< [ \ I W1*U1*V1-C2 

N=0  J=0  K=0 


FIGURE  3-2 

Program  to  Compute  Correlation  Detector  Probabilities 


''TP 
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TABLE  3-1 

False  Alarm  Thresholds 
M *=  WT,  p ■ 0 

x/a^ 


^^l 

2 

5 

10 

20 

50 

100 

.804719 

.599316 

.393138 

.282004 

.200966 

.127726 

.090467 

1.956012 

1.297955 

.775559 

.535128 

.373254 

.234027 

.164992 
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PROBABILITY  OF  DETECTION  (r) 


INPUT  SIGNAL-TO-NOISE  RATIO  (h^  - h,  « h2  ),  dB 

FIG.  3-5  RECEIVER  OPERATING  CHARACTERISTICS,  CORRELATION  DETECTOR, 
FOR  PFA  • 0.01  AND  WT  VARIED  (DETERMINISTIC  SIGNAL) 
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0.9999 


0.999 


0.99 


0.9S 


0.9 


0.8 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 


0.05 

0.02 


■5  0 5 10 


SIGNAL-TO-NOISE  RATIO  (h^  - h22).  dB 

FIG.  3-6  RECEIVER  OPERATING  CHARACTERISTICS,  CORRELATION  DETECTOR, 
FOR  PFA  - 0.0001  AND  WT  VARIED  (DETERMINISTIC  SIGNAL) 
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time  Is  200  seconds.  Solution:  Prom  Figure  3-5»  the  vrr-20  curve. 


Iwe  find  that  y".95  corresponds  to  h “-l.l  dB.  Referred  to  one  hertz 

[bandwidth,  this  value  gives  a detection  threshold  of  -11.1  dB. 

EXAMPLE  3-2.  Problem:  For  a correlation  detector  with  VTT-IO,  what 

-4 

Pjj  can  be  achieved  by  a SNR  of  zero  dB  If  a*10  ? Solution:  Prom 

Figure  3-6,  Pp».56. 

example  3-3.  Problem:  What  Integration  time  Is  required  to  yield 

-4  2 

Pp».9  when  a-10  , h *-3  dB,  and  W»2Hz?  Solution:  Prom  Figure  3-6, 

2 

we  see  that  the  point  (Y»h  ) ■ (.9»-3)  falls  between  the  WT*50  and 
WT-100  curves,  say  at  WT*60.  Therefore,  T should  be  set  to  about 
3.0  seconds. 

To  compute  ROC  for  the  correlation  detector  when  the  signals 

are  Gaussian,  the  procedure  of  Figure  3-4  Is  used,  except  that  with 

h -h  "0,  the  summation  is  finite  (Ll»L2«»0  In  program).  The  signal 
3 *♦ 

parameters  are  Included  by  using  the  notation  (3-17)  In  the  program, 
with  no  noise  correlation  (p-0)  and  perfect  signal  correlation  (r**!). 


p In  program 


h h 
1 2 

2 

1 +h 


-2  In  program 


2 2 2 
a in  program  -►  o (1+h  ) In  program. 


Figures  3-7  and  3-8  give  ROC  for  a Gaussian  signal  case,  assuming 

-2  -4 

h -h  “h,  when  a*10  and  10 
l 2 

EXAMPLE  3-4.  Problem:  What  Increase  In  Input  SNR  Is  required  to 

maintain  Pjj“.95  for  a correlation  detector  with  Gaussian  signals 

-2  -4 

and  noise.  If  P^^  Is  decreased  from  10  to  10  , and  WT-2? 


(3-23) 


1 
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-10  0 10  20  30 

INPUT  SIGNAL-TO-NOISE  RATIO  jjg 


FIG.  3-7  RECEIVER  OPERATING  CHARACTERISTICS,  CORRELATION  DETECTOR, 
FOR  PFA  > 0.01  AND  WT  VARIED  (RANDOM  SIGNAL) 
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FIG.  3-8.  RECEIVER  OPERATING  CHARACTERISTICS,  CORRELATION  DETECTOR, 
FOR  PFA  > 0.0001  AND  WT  VARIED  (RANDOM  SIGNAL). 
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Solution:  Prom  Figure  3-7 » ■ 9.8  dB  for  WT  - 2 and  y ■ .95,  In 

Figure  3-8,  for  the  same  WT  and  y*  - 12.3  dB.  Thus  an  Increase 
of  2.5  dB  Is  required. 

EXAMPLE  3-5.  Problem:  How  does  the  detection  performance  of  a 
correlator  compare  with  that  of  a "sum  and  square"  detector  with 
two  Gaussian  Inputs?  Solution:  Assuming  that  the  signals  are 
perfectly  correlated  (r*l),  the  noises  are  Independent  (p=0),  and  that 
Che  SNR's  In  the  channels  are  equal,  the  SNR  at  the  Input  to  the 
square-law  device  Is  twice  that  In  the  Input  channels.  Thus  we  may 
use  Figures  2-4  and  2-5  If  we  subtract  3 dB  from  the  SNR  Indicated. 
Table  3-3  summarizes  the  comparison  to  be  made,  using  data  from 
Figures  2-4,  2-5,  3-7,  and  3-8.  Prom  the  comparison,  we  observe 
that  the  square-law  Is  better  (requires  a smaller  SNR)  than  the 
correlation  detector  for  the  case  modeled,  but  only  by  about  1 dB  or 
less,  typically. 

I 

i 

i 

I 

I 

I 
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CHAPTER  4 

INCOHERENT  POWER  DETECTOR 

4.1  Detector  Configuration 

f 

i 

I 

i A third  class  of  quadratic  detector  which  can  be  analyzed 

conveniently  is  that  in  which  the  inputs  are  squared  and  then  summed, 
as  illustrated  in  Figure  4-1.  The  output  of  the  lowpass  filter 
is,  in  the  lowpass  detector  bandwidth  case, 

L 

y(t)  = I x2(t),  (4-1) 

i=l  ^ 

where  the  Inputs  {x^}  are  correlated  and  have  different  variances 

in  general.  If  they  were  independent  and  of  equal  variance, 

the  sum  would  be  distributed  x^(L)  and  the  detector  output,  x^(LM)* 
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Rather  than  treat  this  specialized  case,  we  shall  maintain  the  genera- 
lity and  restrict  the  number  of  Inputs  to  L=2,  the  cases  of  larger 
dimension  being  encompassed  by  the  approach  given  in  Part  II. 


For  L=2,  then 


y = 


2 2 

X + X 
1 2 


2 2 

= 2(u  +u  ) , 
3 4 


(4-2) 


where  u and  u are  the  sum  and  difference  terms  defined  before  in 
3 4 


(3*’6),  with  the  moments  given  in  (3-7)»  (3-8),  and  (3-9)»  The 
detector  output  is 


M 


2 . 2 . . 2 


" = ft 


- z + z 
1 2 


(4-3) 


For  both  lowpass  and  narrowband  detector  cases. 


■WTz 


is  x''^(2WT,2WTh^) 
3 


DS 


(4^) 


is  x’^(2WT,2WTh^), 
4 
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so  that  z Is  the  weighted  sum  of  correlated  noncentral  chl-squared 
random  variables. 

4.2  Receiver  Operating  Characteristics 

The  pdf  for  the  output  of  the  square  and  sum  detector  with  two 
Inputs  Is  developed  In  Appendix  C.  Prom  (C-l4),  the  probability 
of  detection  Is 


D 


exp 


(1+p) 


2 

+ h +h 
3 


DS 


X 


m CO 


n=0  m=0 


t 2 

I [WTh 
K=0  *♦ 


n 2 m 
] (WTh  ) 

3 


r_2,\K  (n+WT)j^ 

1-pJ  nlmlKf  ‘^m+n+K+2WT-l 


where  the  equal  input  noise  variance  assumption  has  been  made.  The 
corresponding  false  alarm  probability  is 


- j < p < 1 . 

(^-6) 
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For  deterministic  signals,  the  computer  program  of  Figure  4-2 

was  written  to  compute  (4-5)  when  the  input  SNR's  are  equal  and  the 

signals  are  in  please.  The  program  was  used  according  to  a strategy 

very  analogous  to  that  diagramir:ed  in  Figures  3-3  and  3-^*  False 

— ? -4 

alarm  thresholds  for  a=10  and  10  are  given  in  Table  4-1,  and 

ROC  for  the  same  a are  shown  in  Figures  4-3  and  4-4  for  deterministic 

signals. 


EXAMPLE  4-1.  Problem.:  How  do  the  three  types  of  detectors  we  have 

-4 


considered  compare  as  to  minimum  detectable  signal  (Pp=.5,  o=10  ) ' 

as  WT  varies,  for  two  channels  of  uncorrelated  noise  and  a deterministic 
signal  common  to  the  charnels?  Solution:  Using  Figures  2-7,  3-6, 
and  4-4,  the  comparison  in  Figure  4-5  was  constructed,  showing  that 
in  this  instance  the  coherent  power  or  square-law  detector  is 
consistently  "best"  and  the  incoherent  is  consistently  "worst", 
although  the  difference  in  MDS  is  at  most  1.6  dB. 


For  the  Gaussian  signal  case,  (4-6)  was  computed  using 


2 2 2 2 2 
p=h  /(1+h  ) and  replacing  a with  o (1+h  ) in  the  program,  as  done 


previously  for  the  correlation  detector.  Figures  4-6  and  4-7 


provide  the  ROC  for  the  two-input  sum  and  square  detector  when  a=10 

,-4 


-2 


and  10  respectively. 

EXAMPLE  4-2.  Problem:  Perform  the  same  comparison  as  in  Example 


4-1  for  Gaussian  signals.  Solution:  Using  data  from  Figures  2-5, 
3-8,  and  4-7,  the  comparison  plotted  in  Figure  4-8  shows  much  tbe 
same  result  as  in  the  deterministic  signal  case. 


53 


I 


NSWC/WOL/TR  76-1^8 


» 

t 


( 


I 

I 

I 

I 


i ! 


80  DO=2«M»H/(i+RO) 

82  AO='-2*RO/(l-RO) 

85  IP  D0>0  THEN  100 

86  L1»0 

100  y=X»M/(l+RO) 

105  S=0 
110  Dl=l 

120  FOR  1=0  TO  LI 
130  Al=BO=Bl=CO=l 
140  IF  1=0  THEN  160 
150  D1=D1»D0/I 
160  FOR  J=1  TO  I+2»M-1 
170  BO=BO»Y/J 
180  B1=B1+B0 
190  NEXTJ 

194  S1=B1 

195  IP  RO>0  THEN  200 
197  L2=l 

200  FOR  K=1  TO  L2 

210  A1=A1»A0/K  - 

230  BO=BO*Y/(K+2»M-l+I) 

240  B1=B1+B0 
245  CO=CO»(M+K-l) 

250  S1=S1+A1»C0*B1 
260  NEXT  K 
265  S=S+S1»D1 
270  NEXT  I 

275  S=S»EXP(-Y-DO)»((l+RO)/(l-RO))  +M 
H = h2 
RO  = p 


X = T/a2 
M = WT 

FIGURE  4-2 


I 

f 


Program  to  Compute  Square-And-Sum  Detector  Probabilities 


r 


54 


NSWC/WOL/TR  76-148 


TABLE  4-1 

False  Alarm  Thresholds,  Square  and  Sum  Detector 


Threshold, 


0 = 10"2 

lO”** 

CM 

1 

0 

rH 

6.638352 

11.756371 

3.2798 

5.02256 

7.95691 

3.0226 

3.75662 

5.23860 

2.7775 

3.18454 

4.10311 

2.6487 

2.80822 

3.39456 

2.5558 

2.49445 

2.8306O 

2.4723 

2.34362 

2.56828 

2.4298 

— 

— 

2.3264 

Normalized  threshold  (^^a^-2 


10~‘* 


6.8988 

5.9569 


OF  DETECTION  (r) 
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-5  0 5 10 


INPUT  SIGNAL-TO-NOISE  RATIO  (h^  = - hf  ),  dB 

FIG.  4-3.  RECEIVER  OPERATING  CHARACTERISTICS,  INCOHERENT  DETECTOR, 
FOR  PFA  - 0.01  AND  WT  VARIED  (DETERMINISTIC  SIGNAL) 


PROBABILITY  OF  DETECTION  ( T) 


-5  0 5 10 


INPUT  SIGNAL-TO-NOISE  RATIO  {h^  = = h|),  dB 

FIG.  A-A.  RECEIVER  OPERATING  CHARACTERISTICS,  INCOHERENT  DETECTOR. 
FOR  PFA  ° 0.0001  AND  WT  VARIED  (DETERMINISTIC  SIGNAL) 
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W 


SYSTEMS  COMPARISON,  DETERMINISTIC  SIGNAL 


FIG.  4-6  RECEIVER  OPERATING  CHARACTERISTICS,  INCOHERENT  DETECTOR, 
FOR  PFA  0.01  AND  WT  (RANDOM  SIGNAL) 
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PROBABILITY  OF  DETECTION  (r) 
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FIG.  4-7  RECEIVER  OPERATING  CHARACTERISTICS,  INCOHERENT  DETECTOR, 
FOR  PFA  > 0.0001  AND  WT  VARIED  (RANDOM  SIGNAL) 
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CHAPTER  5 

APPROXIMATION  TO  THE  DISTRIBUTION  OF  QUADRATIC  FORMS 


Having  treated  In  the  previous  three  chapters  selected  cases  of 
quadratic  detectors,  we  now  turn  to  the  more  general  case  In  which 
the  detector  output  Is  an  arbitrary  quadratic  form  of  the  Input 
samples.  Whereas  before  the  assumption  that  the  samples  are 
Independent  permitted  direct  calculation  of  probabilities,  the  removal 
of  that  assumption  now  will  dictate  resorting  to  an  approximation 
technique.  However,  the  approximation  used  Involves  an  asymptotic 
series.  In  which  any  desired  accuracy  can  be  obtained  by  controlling 
the  number  of  terms. 


In  this  chapter,  the  theory  behind  the  approximation  technique 
known  as  Edgeworth's  series  Is  developed  for  quadratic  detectors. 
Specific  computational  procedures  and  examples  of  receiver  operating 
characteristics  obtained  by  this  technique  are  given  In  Chapter  6. 


5.1  Characteristic  Function 

Let  M samples  of  the  L detector  Inputs  be  assembled  to  form 
a vector: 

C’  - . . .Xj^^;x22»  • • • . ;Xt„,X5„,  . . . ,x^„]  (5-1) 


‘IM*  2M' 


where 
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E{U'}  = 


Ll  I 


M-2 


I2  • • • • 


• • • • 


'M-3 


'M-1 


'M-2 


=I 


(LMxLM) 


(5-3) 


where  each  sub-matrix  is  LxL  and 


= I |E{x^(t)Xj  (t+KAt)}|  I 

= ||R^j(KAt)||  , (5-4) 

Rij(t)  being  the  cross-correlation  function  between  input  channels 
i and  J . If  a deterministic  signal  vector  p is  defined  analogous 
to  C in  (5"1),  the  Joint  probability  density  function  for  the 
components  of  C is  given  by 

P(C)  = [(2Tr)^”detj;]”^  exp{  -|(C-y)’^  ^(5-m)}  (5-5) 


in  the  lowpass  detector  case,  and  by 
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Adopting  the  detector  model  of  Figure  1-1,  we  postulate  a 
detector  output  of  the  form 


t 


, w 

z ■ ^ I Q5  , (Q  LXL) 


1-1  1 1 


with 


^1  “ ^^11**21* • * * **L1^ 


Alternately,  this  output  can  be  written 


^ - E'8„5  or  + ^;e^E3 


(lowpass) 


(narrowband) 


in  which 


°L'  * 


^L-  • 


(5-7) 


(5-  8) 


(5-9) 


(5-10) 


an  LMxLM  block  diagonal  matrix.  (0^^  denotes  an  LxL  null  matrix.) 

The  characteristic  function  for  the  distribution  of  the  detector 
output  z is  given  by 


♦^(ip) 


E{lpz}  = EdpC’Qjjj?}  (lowpass) 
^^¥^cVc^^^^^sVs^  (narrowband) 


(5-11) 
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Using  the  results  of  Middleton  [18,  Chapter  17],  who  likes  to  refer 
to  this  distribution  "generalized  x* we  have 

t^Clp)  * (det  A)“^/^txp{- (5-12) 
In  the  lowpass  case,  and 


*2 (ip)  = (det  A)"*exp{- )p^. 


In  the  narrowband  case,  using 


Ilj/j-lP^Qm’  narrowband 


(5-13) 


(5-14)  I 


Now,  although  these  expressions  are  exact,  the  corresponding  probability 
density  functions 


^ /'^p  4^(lp)e-^P“ 

* gp 


(5-15) 


are  "Intractable"  except  for  LM  = 1,  2 and  the  degenerate  case  where 
according  to  [18].  Thus  another  approach  Is  taken. 
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5.2  Asymptotic  expansion  for  the  pdf. 
We  use  Mlddletonb  formula. 


- A"'  ■ - I,  (I9m> 

(lowpass) 


(5-16) 


and  the  relation 


det  A - exp  {-J  tr[(jQj^)”^]}  (5-17) 


m*l 


(lowpass) 


where 


tr[C]  = X c 


11 


(5-18) 


Is  the  trace  of  the  matrix  C.  Thus  we  can  write  the  characteristic 
function  as 


(^^(Ip)  = exp  aj^dp)*"} 
m*l 


(5-19) 


with 


for  a lowpass  detector,  and 

for  a narrowband  detector.  (5-21) 

For  small  correlation  between  time  samples  (i|)^<<l)  and  a lowpass 
sampling  rate  twice  that  of  the  corresponding  narrowband  detector.  It 
Is  shown  In  Appendix  D that  (5-20)  and  (5-21)  are  equivalent. 
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Putting  together  (5-15)  and  5-19),  we 


have  the  asymptotic  expansion  with  respect  to  the  number  of  variables  LM 


OD 

P^^z)  “ ^ ®*P  ^"^P^  I 

— • m*l 


CD 


dp  exp{-lp(z  - a ) - p^a  } exp  a (ip)"*} 


= — f dp  exp{-ip  ( Ij-  p2/2}exp{Jb  (ip)®}  (5-22) 

2ir/2ir  J \/2ir/  ® 


where 


b„  * a /(2a 

ni  in  2 


(5-23) 


Now,  the  first  term  in  the  Integrand  by  itself  would  yield 


2iT/2a 


^ J dp  exp{-lp^^^\  - 

n — OD  ' 


p2/2}  = (2/ira  )“^  exp{ 
2 


*(^)> 


(5-24) 


in  which  (^(x)  is  the  pdf  of  a Gaussian  variable  with  zero  mean  and  unit 
variance.  The  second  factor  in  the  integrand  of  (5-22)  can  be  recognized 
as  a power  series: 


exp( J b^dp)")  = I i [ I b„(ip)»]' 

m*3  K*(J  m=3 


" * \Kl 

m“3  m“3  n“3 


+ . . . 
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1 + I (ip)"* 


(5-25) 


Since  the  integral  (5-22)  is  actually  a Fourier  transform,  we  recall 


that 


P"i{(lp)^  F[g(z)]}  - (-I)*'  g^^^  (z).  (5-26) 


Thus  we  find  that 


Pz(z)  “ 


O ' o 


/Jir" 

2 


(5-27) 


Hem(*)  being  the  Hermite  polynomial. 

Accurate  representation  of  the  density  function  by  means  of  the 
asymptotic  expansion  (5-27)  requires  that  terms  be  arranged  in  order 
of  magnitude  with  respect  to  dimension  LM.  Truncation  of  the 
series  then  yields  an  error  no  larger  than  the  absolute  value  of 
the  first  discarded  term  [16,  S 0.33]. 

A truncation  of  (5-27)  accurate  to  0[  6.5]  is  the  finite  series, 

P»(z)  = — ^ ■ ♦ (z*)  I 1 + b He  (z')  + [b  He  (z')  + b^He  (z*)] 

Z /2r"  ( 3 3 4 4 3 6 


+ [b  He  (z')  + b b He  (z')  + ^ b^ He  (z')] 

55  3U7 


3 4 7 


(5-28) 


+ [b  He  (z')  + (JsbZ  + b b ) He  (z')  + *5  b^b  He  (z') 
66  4 3 5 8 3410 
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where  bracketed  terras  are  of  the  sarae  order  of  raagnltude,  and 


z'  = (z/-  a )/i/2  a . 

1 2 


(5-29) 


5.3  Asyraptotic  expansion  for  the  probability  Integral. 

The  probability  Integral  (or  corapleraentary  curaulatlve 
distribution  function)  for  the  output  of  the  quadratic  detector  is 

P(t)  ® Pr{z>T} 


dz  Pz(z) 


P(t)  = 


/dx  P_ (x/2a  + a ) 

z 2 1 


/TiT 

2 


= Q(t')  + 4.(T’){b  He  (t')  + [b  He  (t')  + h b^ He  (x*)] 

3 2 4 3 3 5 

+ [b  He  (t*)  + b b He  (x’)  + J b^He  (x’)] 

54  346  038 

+ [b  He  (x')  + w + b b ) He  (x’)  + b^  b He  (x’) 
65  4357  34  9 


(5-30) 


again  for  convenience  using 


t'  = (x  - a )//2a  . 

1 2 


(5-31) 
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CHAPTER  6 

COMPUTING  R.O.C.  FOR  GENERAL  QUADRATIC  FORMS 


Although  sufficient  Information  is  provided  in  the  previous 
chapter  to  compute  receiver  operating  characteristics  for  the 
more  general  cases  of  quadratic  detectors,  it  is  appropriate  to 
document  some  of  computational  techniques  which  have  been  used 
with  success.  In  this  chapter  these  techniques  are  described 
and  illustrated  with  example  calculations,  and  results  similar 
to  those  in  Chapters  2-4  are  given  also. 

6.1  Generation  of  Expansion  Coefficients 

A general  outline  of  the  procedures  followed  in  this  chapter 
is  diagrammed  in  Figure  6-1.  As  Indicated  in  the  figure,  the 
first  task  is  compute  the  expansion  coefficients  used  in  (5-28) 
and  (5-30),  based  on  signal,  noise,  and  detector  parameters. 

The  remainder  of  the  procedures  are  useful  for  computing  any 
distribution  for  which  these  coefficients  have  been  generated. 

From  (5-21)  and  (5-23),  expansion  coefficients  for  narrow- 
band  detectors  are  given  by 
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FIG.  6-1.  OUTLINE  OF  COMPUTATION  PROCEDURES. 
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b = a /(2a  m > 3, 

m m 2 “ 


(6-2) 


where  I and  Qj^j  take  the  form  shown  in  (5-3)  and  (5-10),  respectively , 

To  begin,  the  user  chooses  values  for  L,  M,  J,  and  Q. 

The  L X L matrix  Q is  chosen  to  reflect  the  detector  configuration. 
For  example,  for  a (two-input)  correlator. 


[0  .. 
Lk  OJ  . 


(6-3) 


Construction  of  the  LM  x lm  covariance  matrix  involves  assuming 

the  form  of  the  submatrices  J . Here  we  shall  assume  "pink  noise," 

k 


RlJ(kAt)  = Uhl* 


(6-4) 


where  correlation  between  time  samples  decreases  geometrically 
with  the  time  separation  of  the  samples.  With  this  given,  the 
sub-matrices  take  the  form 


Ik  ■ *"lo. 


(6-5) 


which  makes  it  rather  simple  to  generate  as  performed  by  the 
program  in  Figure  6-2,  which  computes  ][,  and  ^Qj^,  given 
L»  M,  Q,  and  ip.  Figure  6-3  illustrates  the  output  of  the 
program  for  L = M = 2,  Q as  in  (6-3),  ^ = .1,  and 


^0  = t ?]  . 

Once  the  signal  components  and  p_  are  specified, 

c s 

computation  of  the  expansion  coefficients  is  straightforward. 
For  the  stationary  signal  case,  which  was  assumed  throughout 


(6-6) 
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5 

BASEl 

10 

DIM  Q(40,40),Z(40,40) 

,V(10,10),W(10,10) 

15 

DIM  R(40,40) 

20 

M-2 

30 

L-2 

40 

W(l,l)-W(2,2)-1  , r 
W(l,2)-W(2,l)-.5^ 

60 

80 

P-.l  ° 

90 

V(l,l)-V(2,2)-0  , o 

V(l,2)-V(2,l)-.5  ^ ^ 

100 

105 

MAT  Z-ZER(L*M,L»M) 

106 

MAT  Q»ZER(L»M,L»M) 

107 

MAT  R-ZER(L»M,L»M) 

no 

FOR  I-O  TO  M-1 

120 

FOR  J-1  TO  L 

130 

FOR  K-1  TO  L 

140 

J1-J+I»L 

150 

K1-K+I»L 

160 

Z(J1,K1)-W(J,K) 

170 

Q(J1,K1)-V(J,K)/M 

180 

FOR  Il-l  TO  M-l-I 

190 

K2-K1+I1»L 

200 

Z(J1,K2)-(P+I1)»Z(J,K) 

210 

Z(K2,J1)-Z(J1,K2) 

220 

NEXT  11 

230 

NEXT  K 

240 

NEXT  J 

250 

NEXT  I 

260 

PRINT  ” 

MATRIX 

Z" 

270 

PRINT 

280 

MAT  PRINT  Z 

290 

PRINT" 

MATRIX 

Q" 

300 

PRINT 

310 

MAT  PRINT  Q 

320 

MAT  R=Z»Q 

330 

PRINT" 

MATRIX 

Z»Q" 

340 

PRINT 

350 

MAT  PRINT  R 

500 

STOP 

600 

END 

FIGURE  6-2 

Program  to  Construct  Matrices 
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1 

MATRIX  Z 

.5 

.1 

.05  “ 

.5 

1 

.05 

.1 

.1 

.05 

1 

.5 

.06 

.1 

.5 

1 _ 

0 

MATRIX  Q 

.25 

0 

0 “ 

.25 

0 

0 

0 

0 

0 

0 

.25 

0 

0 

.25 

0 J 

.125 

MATRIX  Z*Q 
.25 

.0125 

.025  “ 

.25 

.125 

.025 

.0125 

.0125 

.025 

.125 

.25 

.025 

.0125 

.25 

.125  _ 

[ FIG.  6-3  EXAMPLE  OF  COMPUTED  MATRICES. 

[ 


) 

i 

£ I 


Qm 
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m 


chapters  2-^1, 


“cl  ■ 
“si  • 


1 ■ J + nL 

J - 1,  2,  ...L  ^ 

n « 0,  1,  2,  M-1 


(6-7) 


If,  as  a reference  case,  we  suppose  that 


- e 


all  J, 


(6-8) 


then  (6-1)  takes  the  simple  form 


with 


y'  = (1,  1,  , 1).  (LM  elements)  (6-10) 

The  program  of  Figure  6-4  Implements  the  expressions  (6-9)  and 
(6-2),  and  for  the  simple  matrices  of  Figure  6-3,  gives  the  expansion 
coefficients  listed  In  Table  6-1. 

"example  6-1.  Problem:  Find  the  expansion  coefficients  for  a sum- 
and-square  detector  with  four  Inputs  and  WT  = 10,  assuming  no 
signal,  an  Intersample  correlation  factor  of  i|»  * .12,  and  an 
Input  covariance  matrix 


(6-11) 


■ 1 

.24 

.05 

.01 

.24 

1.5 

.3 

.05 

.05 

.3 

1.5 

.24 

.01 

.05 

.24 

1 

Solution:  For  this  configuration, 

n 1 


1 1 

1 1 

1 1 
1 1 


(6-12) 


r 
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BASEl 

DIM  0(^0,40),Z(40,40),V(10,10),W(10,10) 

DIM  R(40,40),U(40,1),Y(40,1),F(1,40),A(12) 

DIM  1(1,1) 

M-2 

L-2 

L1-L»M 

W(1,1)=W(2,2)-1 

W(l,2)»W(2,l)-.5 

P*.l 

S»1 

V(l,l)=V(2,2)-0 
V(1,2)*V(2,1)».5 
MAT  Z»ZER(L»M,L»M) 

MAT  Q»ZER(L»M,L»M) 

MAT  R=ZER(L»M,L»M) 

FOR  I-O  TO  M-1 
FOR  J»1  TO  L 
FOR  K-1  TO  L 
J1=J+I»L 
K1=K+I»L 
Z(J1,K1)=W(J,K) 

Q(J1,K1)=V(J,K)/M 
FOR  11=1  TO  M-l-I 
K2=K1+I1»L 

Z(J1,K2)*(P+I1)»Z(J,K) 

Z(K2,J1)=Z(J1,K2) 

NEXT  II 

NEXT  K 

NEXT  J 

NEXT  I 

MAT  R«Z»Q 

MAT  U-CON(LI,I) 

MAT  y-CON(LI,I  405  A(K)=A(K)/2 

MAT  Y-Q»U  410  IF<K  3 THEN  430 

MAT  F-C0V(1,L1)  420  A(K)=A(K)/(2»A(2) )+(K/2) 

MAT  F-TRN(Y)  430  MAT  Q=(1)»Z 

MAT  Z-(1)»R  440  MAT  Z=Q»R 

MAT  Q=IDN(L1,L1)  450  PRINT  K,A(K) 

FOR  K=1  TO  12  460  NEXT  K 

MAT  Y»Q»U  500  STOP 

MAT  T=P*Y  600  END 

A(K)-0 

FOR  J»1  to  LI 
A(K)-A(K)+Z(J,J) 

NEXT  J 

A(K)-S»I(1,1)+2»A(K)/K 

FIGURE  6-4 

Program  to  Compute  Narrowband  Coefficients 
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TABLE  6-1 


EXAMPLE  NARROWBAND  CORRELATOR  COEFFICIENTS 

o2=o2»l,  p».5,  S »S  =1,  0=0  ; L*M»2 

12  12  12 


m 


or 

m m 


1 

2 

3 

!♦ 

5 

6 

7 

8 

9 

10 

11 

12 


.36^062  ^ ^ 
.193056 
8.62109E-2 
3.91863E-2 
I.82383E-2 
8.5753OE-3 
^.06278e-3 
1.93408E-3 
9.2387IE-4 
^♦.42390E-4 
2.12217E-4 


TABLE  6-2 


EXAMPLE  SQUARE-LAW  COEFFICIENTS 

(6-11);  t|«=.12. 


m 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


COVARIANCE  MATRIX  GIVEN  BY 
L=4,  M=10 


6.78  , „ 

2.35877^  ^ 2 

.10948 

2.75998E-2 

7.57808E-3 

2.20737E-3 

6.71795E-4 

2.11502E-4 

6.83994E-5 

2.26049E-5 

7.60434E-6 

2.59604E-6 
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and  in  the  program  of  Figure  6-3  is  defined  by  the  statements 

90  MAT  V = CON(4,4) 

100  (deleted). 

The  covariance  matrix  is  entered  by  the  statements 

40  MAT  W = ZER(4,4) 

60  MAT  READ  W 
510  DATA  1,  .24,  .05,  .01 
520  DATA  .24,  15,  .3,  .05 
530  DATA  .05,  .3,  1.5,  .24 
540  DATA  .01,  .05,  .24,  1 
Other  changes  to  the  program  required  are 
20  M = 10 
30  L = 4 
80  P = .12 
85  S = 0. 

The  coefficients,  as  computed  by  the  program,  are  given  in  Table  6-2. 

As  a quick  check,  note  that  the  mean  value  of  the  detector  output, 

given  by  a , agrees  with  the  correct  value,  which  in  this  case  can 

easily  be  calculated  by  summing  all  the  elements  of  Finally, 

M 

since  \li  <<  1,  these  results  apply  for  both  lowpass  and  narrowband 
detectors . 
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6.2  Computation  of  PDF  and  CDF 

After  obtaining  the  coefficients  for  the  asymptotic  expansion 
(5-19)  for  the  characteristic  function  of  the  detector  output's  i 

probability  distribution,  it  is  straightforward  to  compute  j 

either  the  probability  density  function,  using  (5-28)  or  the 
cumulative  distribution  function  using  (5-30).  Common  to  both 
of  these  expressions  are  the  Kermite  polyrKjmials , which  are 
given  by 

He  (x)  = (-1)"  [ ■ <f)  (x)]/<()(x) 

dx'* 

= xHe^_^(x)  - (n-1)  He^_^(x), 

with  HeQ(x)  “ 1,  He^(x)  = x. 

Continuing  the  example  of  the  correlator  with  given 

by  (6-6)  and  Q by  (6-3)  and  with  h^  = h*  = 1 (S  = /2^  in  (6-8))  and 

1 2 

ij/  = .1,  Figure  6-5  shows  the  detector  output  PDF  for  different 
values  of  M,  the  number  of  samples.  Two  effects  are  evident 
in  the  figure  as  M Increases:  the  variance  of  the  distribution 
decreases  (the  peak  becomes  narrower),  and  the  distribution  begins 
to  approach  symmetry  about  the  mean  value  ( = 1.5),  indicating 
a slow  convergence  to  a Gaussian  distribution.  These  effects 
are  perhaps  more  obvious  in  Figure  6-6,  in  which  the  probability 
Integral  (5-30)  is  plotted,  using  the  following  approximation  [l] 


(6-13) 


for  the  Gaussian  probability  integral: 

Q(x)  ( (|)(x)  [d  y + d y^  + d y*],  x > 0 (6-14) 

J 12  3 — 

“ (l  - Q(-x),  X < 0 


y = 1/(1  + .33267X) 

d = .4361836 
d'*  -.1201676 
d*“  .9372980. 

3 
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CORRELATOR  OUTPUT  PDF,  p(z) 


CORRELATOR  OUTPUT  (z) 


FIG.  6-5.  CORRELATOR  PROBABILITY  DENSITY  FUNCTION. 
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EXAMPLE  6-2.  Problem:  What  effect  do  inter-sample  and  Inter-channel 
correlation  have  on  the  distribution  of  an  incoherent  power 
detector  with  five  inputs  and  M « 6?  Solution;  In  terms  of  the 

notation  we  have  adopted,  L ■ 5 and  Q ■ I , The  effect  of  the 

s 

inter-sample  correlation  is  found  by  varying  the  parameter  4). 

Inter-channel  correlation  is  modeled  by  the  covariance  matrix 

Let  us  assume  that  the  sensors  feeding  the  five  channels  are 

so  oriented,  and  the  noise  field  is  such  that  the  elements 

of  y are  given  by 

° (1,  i - J 

E(n.n.)  = < (6-16) 

^ ^ ( P,  i i 

With  these  parameters  specified,  the  computation  of  the  pdf 
and  probability  Integral  is  accomplished  by  modifying  the  program 
of  Figure  6-3  and  by  the  addition  of  the  statements  listed  in 
Figure  6-7.  The  probability  Integral  for  this  case  is  shown 
in  Figures  6-8  and  6-9  for  ij'  varied  and  p varied,  respectively. 

It  is  shown  in  these  figures  that  positive  correlation  increases 
the  dispersion  (variance)  of  the  distribution  and  that  negative 
inter-sample  correlation(il^<0)  can  decrease  the  variance  for  small 
values,  but  has  a degrading  effect  in  the  sense  that  convergence 
to  a Gaussian  distribution  is  retarded  (as  evidenced  by  an  increase 
in  curvature  in  Figure  6-8  for  4'<0)» 

6.3  Computation  of  ROC 

Given  the  ability  to  compute  the  probability  integral 
for  the  types  of  quadratic  detectors  we  have  been  considering, 
it  is  relatively  simple  to  obtain  receiver  operating  characteristics 
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440  MAT  Z=Q«R 
450  NEXT  K 

460  REM  COMPUTE  PPF,  PROB  INTEGRAL 
465  PRINT  "X","P(X)","Q(X)" 

470  X1=0 
480  X2=4 
490  X3=.5 

500  FOR  X=X1  TO  X2  STEP  X3 
510  X4=(X-A(l))/SOR(2»A(2)) 

520  H(2)=X4+2-1 

530  H(1)=X4 

540  FOR  K*3  TO  12 

550  H(K)=X4»H(K-1)-(K-1)»H(K-2) 

560  NEXT  K 
570  GOSUB  1000 
590  GOSUB  1200 
600  PRINT  X,P2,P3 
605  NEXTX 
610  STOP 

1000  REM  PDF  SUBROUTINE 

1010  P2-1+A(3)*H(3)+A(4)»H(4)+A(3)+2»H(6)/2+A(5)*H(5) 

1020  P2=P2+A(3)*A(4)»H(7)+A(3)+3*H(9)/6+A(6)»H(6) 

1030  P2=P2+(A(4)+2/2+A(3)*A(5))*H(8)+A(3)+2*A(4)»H(10)/2 

1040  P2=P2+A(3)+4»H(12)/24 

1050  P2=P2*EXP(fX4+2/2)/(2*SQR(Pl»A(2))) 

1060  RETURN 

1200  REM  PROB  INTEGRAL  SUBROUTINE 

1210  P3=A(3)*H(2)+A(4)»H(3)+A(3)+2»H(5)/2+A(5)*H(4)+A(3)*A(4)«H(6) 

1220  P3=P3+A(3) +3*H(8)/6+A(6)»H(5)+(A(4)+2/2+A(3)*A(5))*H(7) 
1230  P3=P3+A(3) +2»A(4)«H(9)/2+A(3) +4»H(11)/24 
1235  IF  X4<0THEN  1244 
1240  T1=1/(1+.33267*X4) 

1242  GO  TO  1250 

1244  T1=1/(1-.33267*X4) 

1250  P5=.436i836»T1-.1201676»T1+2+.937298*T1+3 
1252  IF  X4>0  THEN  1256 
1254  P5=EXP(X4+2/2)*SQR(2*P1) -P5 
1256  P3=P5+P3 

1260  P3=P3*EXP(-X4f2/2)/SQR(2*Pl) 

1270  RETURN 
2000  END 

FIGURE  6-7 


Program  to  Compute  PDF  and  Probability  Integral 


I 

I 

1 

I 
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r 


I 


for  them.  As  in  chapters  2-M,  the  first  step  is  to  calculate 
false  alarm  thresholds.  The  probability  integral  is  then  evaluated 
at  these  thresholds  for  various  signal-to-nolse  ratios. 

This  report  is  concluded  with  an  illustrative  system 
comparison  problem. 

EXAMPLE  6-3.  Problem:  Comptire  the  detection  performance  of  the 
quadratic  detectors  shown  in  Figure  6-10,  when  a Gaussian  signal 
is  common  to  the  four  channels,  and  the  inter-channel  noise  and 
inter-sample  correlations  are  zero. 

Solution:  Under  the  conditions  stated,  L = 4,  M = 10, 

= 0,  and 


1 + h^  h^  h^ 


h^  1+h^  h^  h^ 
h^  h^  l+h^  h^ 


h^  h^  h^  i+h^ 


(6-17) 


For  the  three  detectors,  the  quadratic  forms  they  implement  are 
characterized  by  the  matrices 


0 .5  0 0 

.5  0 0 0 

0 0 0 .5 

0 0 .5  0 

1111 

1111 

1111 

1111 


(Correlator)  (6-l8) 


(square-law)  (6-19) 


f 

5 

1 

j 


1 
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(Incoherent ) 


(6-20) 


After  varying  the  argument  of  the  probability  Integral  with  h*=0 
to  determine  the  false  alarm  thresholds  for  a = 10“*  and  lO”** 
for  the  three  detectors;  These  values  of  argument  were  fixed  while 
h*  was  varied,  producing  the  ROC's  of  Figure  6-11  and  6-12. 

Prom  these  curves  it  is  evident  that,  for  the  signal  and  noise 
correlations  stipulated,  the  square-law  configuration 
requires  about  2 dB  less  for  detection  than  the  correlator,  which 
in  turn  requires  about  .3  dB  less  than  the  Incoherent  power 
detector. 


INPUT  SIGNAL-TO-NOISE  RATIO,  dB 

FIG.  6-11  FOUR-INPUT  DETECTOR  COMPARISON,  PFA  ° 0.01,  RANDOM  SIGNAL  (WT  > 10) 
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FIG.  6-12.  FOUR-INPUT  DETECTOR  COMPARISON,  PFA  = 0.0001 , RANDOM  SIGNAL  (WT  - 10) 
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APPENDIX  A 


APPROXIMATION  TO  THE  NONCENTRAL  x*  DISTRIBUTION 


In  [12]  it  is  shown  that  when  the  noncentral  chi-squared  random 
variable  X,  with  v degrees  of  freedom  and  noncentrality  parameter 


X,  is  transformed  into  the  variable  Y by  the  relation 


(A-1) 


the  variable  Y very  closely  approximates  a Gaussian  variable  with 
unit  variance  and  mean  equal  to 


E(Y)  = /X  + Js  Cv  - 1). 


(A-2) 


Thus,  the  probability  integral  for  X is  approximately 


Q(xlv,X)  = Q(/x  -"17  - A + \{X-  rj),  (A-3) 

where  Q(x)  is  the  Gaussian  probability  integral.  Defining  a 

normalized  detection  threshold  d^  analogous  to  the  false  alarm 

threshold  d , we  have 
a’ 

d = - ;5(2WT  - 1)  - /2WTh2  + i5(2WT  - 1)  (A-4) 

Y i“a 


Prom  (2-13), 


s 2/WTd  + 2WT,  and 
^i-a  a * 


(A-5) 


d dHl  d (7d~  ; 

h^(Y,  a,  WT)  = — + W — 2L  + 1 + (A-6) 

/Vf  2WT  ^ WT  2WT 


A-1 
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The  accuracy  of  this-  approximation  to  the  noncentral  chl- 
squared  distribution  was  checked  against  the  table  of  [6] 
for  a » lO"*,  with  the  result 


X Y yapprox.  ^dev. 


2 

8.190 

.5 

.4986 

-.28 

13.881 

.8 

.7998 

-.02 

27.415 

.99 

.9902 

+ .02 

4 

10.232 

.5 

.4974 

-.52 

16.749 

.8 

.7997 

-.04 

31.794 

.99 

.9904 

+ .04 

10 

14.126 

.5 

.4961 

-.77 

22.177 

.8 

.7995 

-.07 

40.021 

.99  .. 

.9905 

+ .05 

20 

18.451 

.5 

.4957 

-.87 

28.162 

.8 

.7994 

-.08 

49.027 

.99 

.9907 

+ .07 

40 

24.542 

.5 

.4955 

-.90 

36.550 

.8 

.7993 

-.08 

61.572 

.99 

.9908 

+ . 08 

100 

36.614 

.5 

.4958 

-.83 

53.103 

.8 

.7994 

-.08 

86.175 

.99 

.9908 

+ . 08 
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APPENDIX  B 

DERIVATION  OP  CORRELATOR  PDF 


As  indicated  in  the  text,  the  output  of  the  correlation  detector 
with  deterministic  signal  input  is  the  difference  of  two  scaled 
noncentral  chi-squared  variables: 


where,  using  M » WT 

2Mz 

V 


z - z 

I 2 


is  x'*(2M,  2Mh*) 

3 


(B-1) 


(B-2) 


2Mz 

— r*-  is  x'*(2M,  2MhM. 


Under  the  simplifying  assumption  of  equal  input  noise  power 
= 0,  z^  and  z^  are  independent,  and 


4Mz 

0^ 


2Mz  2Mz 

- (1  + p)  - P) 

3 4 

= (1  + p)v  - (1  - p)v  . 

1 2 


(B-3) 


The  characteristic  function  for  the  distribution  of  a non- 
central chi-squared  variable  V with  2v  degrees  of  freedom  and 
noncentrality  parameter  X is,  from  [ll]. 


B-1 
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.-\/l"  (X/2)"  1 

* -L'^rr^  rr^ 


21p)v+n 


Therefore,  the  characteristic  function  for  the  distribution  of 
V ■ 4Mz/o* 

5 (ip)  - E{e^P<l  * ‘>>\)  E (e-^P'l  - P>\)  „ 

^ • OQ  (Mh*)"^  (Mh*)"* 

■ J-o 


P(1  + p)JM+n  [1  + 2ip(l  - p) 
Making  use  of  an  expression  in  [14],  we  may  write 

t \n  p*i  ■ ^\t 


[1  - 2ip(l  + p)]"”"^[l  + 21p(l  - p)]^ 

■x  fsrr'i  W"  (¥)"•■■’ 


^r-2lp 


m+M 

+ I 

k-1 


: r::rif¥re¥) 


M+m-k 


LI  + Ji 
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This  expression  allows  a representation  of  the  pdf  of  v In  terms 
of  chl-squared  pdf's: 


P„(v)  » 35  exp{-M(h*  + hj  I i 

V \ ^ / * 'n»0  m= 


:M(l-p)h§/2]"[M(l-p)hg/2; 
n!  ml 


n+m-1 


— _ /2M+m+n-2-jV  ? , v . 

(m*.-!  Ai^ 


(B-8) 


m+M-1  /2M+m+n-2-k\/  , \k  ^ , 

^ (H.m-1  )fe) 


, V < 0, 


where 


PX^(x|2v)  2 


, X > 0. 


2r(v) 


(B-9) 


Alternately,  to  get  the  double  infinite  series  into  a 
more  recognizable  form,  such  as 


jrn^T  ^mn^''^* 


(B-10) 


n*0  m«0 
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we  define  the  polynomials 


■ k*0  \ ® ' k! 


(B-11) 


whose  properties  are  summarized  in  [13].  This  allows  the  pdf  of  the 
detector  output  to  be  written  as 


Pz(z) 


HM 

“V  p 


(^r 


exp{-M(h^  + h^)} 


I ^ [M(l-P)h§/2]^  [M(H-p)hS/2] 

n n!  m! 


n=0  m=0 


(B-12) 


exp{-2Mz/o2(l+p)}  [4Mz/o2(1  - p^)],  z > 0 


exp{2Mz/o2(l-p)}  gU+”IJ;  [-i»Mz/o2(l  - p2)],  z < 0. 


NSWC/WOL/TR  76-li»8 


APPENDIX  C 

DERIVATION  OP  SQUARE-AND-SUM  PDF 


As  shovm  In  the  text,  the  output  of  the  Incoherent  or  square- 
and-sum  detector  with  two  deterministic  signal  Inputs  can  be 
considered  the  sum  of  two  scaled  noncentral  chi-squared  variables: 


z - z + z 
1 2 


(C-1) 


where,  using  M ■ VTT, 


^ is  x’^(2M,  2Mh2) 

3 • ^ 


(C-2) 


is  x’^(2M,  2Mh2) 

0^  4 

4 


Under  the  simplifying  assumption  of  equal  input  noise  power, 

0 ■ 0.  z and  z are  Independent,  and 

2 1 2 


Nz . Mz  J 

(1  + p)  -5-  + (1  - p)  -r 
o o 

I •> 


6 (1  + p)v  -f  (1  - p)v  . 


(C-3) 
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The  Joint  pdf  of  and  is 


P.  (v  , V ) - H e“^^x 
* 1 2 


1 (»'A  V ) I„  . (/A  V ),  V , V > 0 
M-1  11  M-1  22  1 2 


(C-4) 


where 


2Mh2 

3 


2Mh2 

4 


(C-5) 


I„(x)  Is  the  modified  Bessel  function  of  the  first  kind,  of  order  K. 
K 

Consider  the  following  change  of  variables: 


r,  . JL.  + ilzfiiZ 

1 l+A  ^ 9 


Vj  . - 

Then  the  pdf  of  x and  y is 


X » (l+p)v  +(l-p)v  > 0 

1 2 


-2v  -2x 


1+p  I^p 


1 y 1 


p,u.  • >»  p.c  ifp  * 


(1/8)  exp{-H[x^+  + py  - 


lM-1 

:i-pMy+2x]j  2 


* ^M-1  ^’'^i^l+p  * 


y]>  Im-1 
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for  X > 0 amd  - 


i y 1 O' 


since  X = ■ ■ Is  the  variable  of  Interest,  we  obtain  Its 

pdf  by  Integrating  out  y: 


>^(x)  • f dyp^Cx,  y)  - y dyp^Cxi  -y) 

-2x  0 


fj(x)  y**dy  e'py  [y(y  - a)] 


(C-8) 


using  a - 


f^(x) 


F f exp{-3s[X^X^  5^]}  (C-9) 


By  expanding  the  Bessel  functions  Into  their  series  representations. 


we  have 


p (X)  . f (X)  I I [(l-p)Xi/8]"  [(l|p)X,/8r 

* * n»0  m*0  nl (n+M-1) ! ml (m+M-1) 1 


y “ ay  y’’*"-^(p  - y)"*"'!. 


(C-10) 


with 


(0-11) 


! 
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Prom  [l6],  formula  3 •383.1,  we  learn  that 


y“  dy  y”+«-l  („  - y)"+«-l 

0 

_ (m+M-1) ! (n+M-1)  1 ..  m+n+2M-l  ™ /^n 

“ V-fn4-2Mil'yi  “ ^P^[m+Mim+n+2M;-pa/2], 


(C-12) 


where  ) Is  the  confluent  hypergeomefcric  function.  After 

substituting  Into  (C-10),  we  obtain  the  pdf  of  the  detector  output  Z 


- p? 


. {¥) 


? ; [MhUUl 

n-O  1.=0 


■o)]"  (MhS )”  rMz/o^dti,)]"*"' 
m!  (n+m+2M-l) I 

(C-13) 


X P^M+n;  2M+ni+n;  -2pMz/a^(l-p^)] , z ^ 0. 

This  expression  can  be  shown  to  be  a special  case  of  that  given 
In  [17],  page  59* 

By  expanding  the  hypergeometric  function,  the  probability  Integral 
is  found  to  be 

P(t)  -f  p,(z)dZ  exp{-M[jTn7jy  + + hj]) 


(expression  continued) 
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‘Jo  nlo  Jo  t""* 


-P/1 


Special  cases.  (1)  For  » 0, 

3 ^ 


p( 


•>  ■ H”  •■ 


Mt/o^CI+P ) 


®k+2M-l^”^  a^Cl+p)}. 

-1  < p < 4 


or 


P(0  -fe)"  e-M^/p^d-p) 


(2)  For  p ■ 0 (Input  noise  uncorrelated). 


- j < p < 1. 


P(0  . exp(-M(  ^ * hUhdlJ  I 


n-0  in»0 


(C-17) 


C-5 
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APPENDIX  D 

EQUr/ALENCE  OP  LOWPASS  AND  NARROWBAND  DISTRIBUTIONS 


Given  a narrowband  detector  sampling  rate  of  B samples  per 
second  and  a lowpass  detector  sampling  rate  of  2B  samples  per 
second.  For  T seconds,  the  number  of  samples  M In  the  expression 
for  the  narrowband  expansion  coefficients  (5-21)  Is  M ■ BT, 
and  for  the  lowpass  coefficients  M ■ 2BT  In  (5-20). 

For  M ■ 2BT,  from  (5-3)  and  (5-10), 


(2BT) 


and 


^2BT  “ ^ 


Now,  If  we  write 


E 

(BT) 



-e- 

td 

W 

E 

(BT) 

1 

1 

1 

1 

1 

1 

E 

(BT) 

%T 

1 

1 

1 

1 

_i 

^(BT) 

°(BT) 

1 

1 

> 

1 

1 

*^BT 

(D-1) 


(D-2) 


te  Q T -.L 

\ (2BT)  2BTy  2"^ 


JE  Q ) 
'm  (bT)  BT 


m 


(BT)  BT 


m 


, (E  Q )*" 

fm  (BT)  BT 


®m^^(BT?BT^  ^ 


(D-3) 


D-1 
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it  Is  true  that 


(D-4) 


nrp  um 

since  e = 1 and  f = . Thus  e^^^  is  of  the  order  1 + 0(V»  ), 

BT 

and  fjjj  is  0(\|<  ),  with  the  consequence 


_1 

2ni-x 


tr 


X 


[l+0(i|^®^)].  (D-5) 


Therefore,  under  the  assumptions  made  here,  the  noise  power  terms 

BT 

in  (5-20)  and  (5-21)  are  equivalent  for  \()  <<1. 

Equal  signal  power  for  the  two  cases  of  bandwidth  is  represented 
by 

yi(tj)  = Aj^(tj)  cos[wtj  + ^^(tj)]  (narrowband) 

Pl(tj)  = (lowpass) 

If  the  signals  are  stationary  over  the  entire  integration  period, 

lei 


and  also 


NSWC/WOL/TR  76-148 


Q fz  Q 

(2BT)  2BT  V (2BT)  2BT 


(2BT) 


-=^  y'  Q Q 

gin-i  (grpj  BT^^BT)  BT/  (BT) 


X (e  + f ). 

**  ^ m-i  m-i 


(D-8 


Thus,  In  the  sense  of  (D-6)  the  signal  terms  In  (5-20)  and 

BT 

(5-21)  are  equivalent  to  order  0(i<)  ) for  the  assumptions  which 

have  been  made. 
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